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on different targets using a non perturbative time-dependent quantum cal- 
culation. The evolution of the neutron halo wave function shows an emis- 
sion of neutron at large angles for grazing impact parameters and at for- 
ward angles for large impact parameters. The neutron angular distribu- 



te ' tion is deduced for the different targets and compared to experimental data. 

We emphasize the diversity of diffraction mechanisms, in particular we dis- 
cuss the interplay of the nuclear effects such as the towing mode and the 

H 

Coulomb break-up. A good agreement is found with experimental data. 
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I. INTRODUCTION 

A major discovery of the last decade in nuclear physics is the observation of halo nuclei [|l|. 
The presence of these extended systems has been uncovered by break-up studies. However, 
an unambiguous interpretation of these reactions requires a deep understanding of reaction 
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mechanisms such as diffractive effects [[|. In particular, the interplay between nuclear and 
Coulomb dissociation is of major importance. 

On the experimental side, measurements of the neutron angular distribution have been 
performed for the one-neutron break-up of n Be on Au, Ti and Be targets || at 41 MeV 
per nucleon. They clearly present two contributions. One peaked at small angles (below 
10 degrees) which strongly varies with the target ; for the Au target the cross section is as 
high as 50 barns per steradian at small angles whereas for the Be target it is around 1 barn 
per steradian. A second one located at large angles which shows a lesser dependence with 
the target ; the cross section around 30 degrees is about 100 mbarn per steradian for all 
studied targets. The small angle region was understood through models based on Coulomb- 
excitation theory to come from Coulomb dissociation [|3] and the large angle emission was 
thought to be a consequence of nuclear break-up. Perturbative calculations ||,|5| including 
both interactions or non perturbative calculations ||[7j which only account for the Coulomb 
field have been performed in the last years. Only recently the time dependent Schrodinger 
equation was solved || and applied to the break-up of n Be using different numerical tech- 
niques than the one we have used |J. In this later paper, the authors predict the presence 
of neutrons emitted with a high angular momentum in the 11 Be frame which corresponds to 
the neutrons measured at large angles that were mentioned above. This prediction, which 
the authors do not observe when they use the eikonal approximation, is in agreement with 
our calculation as it will be discussed further in this paper. At the same time our paper was 
in submission, we noted the work of ref. [|10j that also presents a method of solving the time 
dependent Schrodinger equation and applies it to the break-up of n Be and 19 C to calculate 
the relative energy between the emitted neutron and the core. 

In our article, we present a non-perturbative model which accounts both for the Coulomb 
and the nuclear effects through the resolution of the full time-dependent Schrodinger equa- 
tion onto a cartesian mesh. Our calculation shows that the emission of neutrons at large 
angles is due to the interaction of the particle with the nuclear potential. Such an emission 
has already been observed in reactions between stable nuclei |lTJ where the nucleon origi- 



nated from the target and not the ejectile as in the case of n Be break-up. This mechanism 
was called "towing mode" as the particle was pulled out from the target and towed by 
the projectile for a short while. Through the resolution of the time-dependent Schrodinger 
equation we reproduced this emission to the continuum with specific angles and energies in 
agreement with the measurements ||. 

In the following, we will show through our calculation that both Coulomb and nuclear 
fields play an important role in the case of the neutron break-up of halo nuclei. 



II. DESCRIPTION AND INPUTS OF THE MODEL 

We use the calculation presented in a previous paper || and apply it to the reactions 
197 Au, 48 Ti, 9 Be ( 11 Be, 10 Be + n). This calculation describes the wave function distortion 
of an initially bound particle as it passes by the potential induced by the reaction partner. 
This is performed in the framework of independent particles. 

The wave function of the neutron halo is deduced from the potential found by N. Vinh 
Mau to reproduce the inversion of the 2s and the lp states in the 11 Be nucleus for which 
a derivative of the Wood-Saxon potential is added at the surface. The potential between 
the neutron and the 10 Be core reads 

V Be (r) = — ^U + 16.a &2 ~ (1) 
1 + e « h + e— 

where the diffuseness a is 0.75 fm and the coefficient a is equal to -10.56 MeV for the 
2s state. The radius R equals to 1.27 x A 1 / 3 fm. Numerically the initial wave-packet is 
obtained by diagonalizing the one-body Hamiltonian in spherical coordinates in a sphere of 
30 fm radius with a space-step of Ar =0.02 fm. The depth Vq was taken to be -40 MeV. The 
wave function calculated in spherical coordinates for n Be is then mapped onto Cartesian 
coordinates in order to calculate the dynamical evolution. Special attention has been given 
to the purity of the ground state 11 Be neutron halo by performing an additional imaginary 
time evolution on the Cartesian network, after which the 2s state is bound with -0.503 MeV. 
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The dynamical evolution of the system is given by the following single particle 
Schrodinger equation which reads in the r space 

ihj t v(r) = (^-A + V T (r - rr(f)) + V Be (r - r fle (t))) <p(r) (2) 

where Vr and Vs e are the time-dependent potentials between the neutron and the target 
and between the neutron and the projectile respectively. r^(t) and re e (t) correspond to the 
target and the projectile positions respectively. The nuclear potential of the Au and Be 
target, Vr(r), is taken to be of Wood-Saxon shape with a diffuseness a = 0.5 fm, a radius 
R = 1.27 A 1 / 3 fm (A being the mass number of the 197 Au or 9 Be) and of depth adjusted to 
obtain the experimental binding energy of the last neutron. For the Au target, Vq = -49.3 
MeV, giving a 3p state bound by about 8 MeV and for the Be target, V = -40. MeV, 
giving a lp state bound by about 1.66 MeV. For the 48 Ti target, we used the Becchetti 



and Greenlees |13| prescriptions to obtain the real part of the nuclear potential between the 
neutron and the Ti, leading to a depth of V = -41.2 MeV, a diffuseness of 0.75 fm and a r 
value of 1.17 fm. 

In order to take into account the target and the 10 Be core displacement, the evolution 
is performed for each time step using the classical Coulomb trajectory for the center of 
mass motions (ry(t), fr(i), rs e (t), TBe{t)) from a distance of -400 fm along the initial 
velocity axis between the projectile and the target and equal to the impact parameter b 
on the perpendicular axis. This explicit treatment of the core recoil is responsible for the 
Coulomb excitations of the neutron around the w Be core. The numerical method used for 
the trajectory calculation is the Runge-Kutta method. 



We have chosen to use the split-operator method |15|] for the time evolution. This is a 
well known way to solve the time dependent Schrodinger equation (Eq|J) on a 3D cartesian 
lattice. It is faster and as accurate as the methods based on Taylor expansion of the evolution 
operator which are also routinely used to solve time dependent problems. Moreover it allows 
to treat large lattice performing a rather long time calculation. 

The time step used is of 1.7 fm/c on a mesh of 64x64x64 fm 3 with a step size of 0.5 fm. We 



have tested the numerical accuracy of the method with respect to the mesh parameters. The 
calculation is performed in the initial projectile frame and a Galilean frame transformation 
is performed to extract the observed quantities. 

III. RESULTS 
A. Density plots 

The result of these calculations (Fig .[ID is presented for the n Be + Au reaction as the 
probability density integrated over the z-axis, p2 S {x,y) = J <p(x,y, z) 2 dz, for impact param- 
eters b = 10, 12, 15, 20 and 40 fm, after the evolution of the 2s wave function initially in 
the 10 Be potential. 

After the evolution, we subtract the contributions of the wave-functions which are bound 
in the 10 Be potential in order to keep only the emitted part of the wave-function. This 
subtraction is performed by projecting out all the bound states \a > of the neutron in the 
10 Be, leading to the emitted wave function \ip >. 

>= \(f) > -^2\a >< a\(f) > (3) 

a 

The bound wave functions in the 10 Be core potential are the 2s state itself bound by -0.503 
MeV, the Is state bound by -23 MeV and the lp state using the same nuclear potential and 
which is then bound by about -11 MeV. Note that we do not consider the contribution of 
the lp wave function which has been experimentally evidenced at -0.183 MeV since it is not 
a bound wave function of our potential. 
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FIG. 1. Density plots of the 2s wave function in 11 Be after scattering on a 197 Au target for 
impact parameters of 10, 12, 15, 20 and 40 fm. The Coulomb trajectory is represented by the 
dashed line for each impact parameter. These plots are displayed in a logarithmic scale, black 
representing the most intense density. The ratio in density is 3 between each grey area. 

Fig.0 shows the same evolutions as Fig.[l| after the subtraction of these components. A 
sizeable fraction of the wave function is then removed around the 10 Be core position. One 
can see that the small impact parameters (b < 20 fm) are responsible for neutrons emitted 
at large angles in the opposite direction of the core compared to the initial direction. For 
larger impact parameters, neutrons are forward focussed. This small angle emission can 
be understood as a Coulomb dissociation in which the 10 Be core has been shaken by the 
Coulomb field of the target. In our calculation, for large impact parameters, we observe 
a displacement of the neutron wave function compared to the core which might be the 



first oscillation of a soft dipole resonance as it has been suggested in ref. ||16|| . Due to 
the weak binding of the last neutron, the halo has separated from the core, transferring 
little momentum to the neutron which is then emitted along the initial trajectory. The 
calculation shows a weak emission (see Fig.^] right) which must, however, be integrated over 
a large impact parameter domain to obtain the total cross section. 

In order to compute cross-sections, we used the following method. For large impact 
parameters the fraction of wave function emitted is small and thus according to the per- 
turbation theory, to gain time and reduce the error inherent to the method, we performed 
the calculation up to the minimum distance of approach and multiplied the result by two. 



For large impact parameters we have tested that this procedure gives the same results as 
the complete calculation demonstrating that the impact parameters are large enough to 
guarantee a first order perturbation. 




an impact parameter of 10, 12, 15, 20 and 40 fm and after subtraction of the bound wave functions. 

On the contrary, for small impact parameters (below 20 fm) the nuclear break-up cross 
section is large and the nuclear potential spreads the wave function on the whole mesh 
(see Fig.|l] and @) we stopped the calculation at a distance of 20 fm after the projectile has 
passed by the target and checked that the extracted values do not change for a slightly 
longer evolution. Since the nuclear break-up is much larger than the Coulomb break-up 
cross section for the small impact parameters, as will be shown in section B, the fraction of 
Coulomb break-up missed there is negligible. 



B. Interplay between nuclear and Coulomb break-up. 

Since the 10 Be core is detected in the experiment, we assumed a minimum impact pa- 
rameter corresponding to a grazing trajectory. We followed the strong absorption model for 
which the minimum parameter is calculated as 1.4(A^ 3 + A 1 ^) fm yielding to b min = 11 fm 
for the Au target, b m i n = 8 fm for the Ti target and b m i n = 6 fm for the Be target. To obtain 
a total cross section, this calculation was performed for impact parameters running from the 
grazing parameter, to an impact parameter where the fraction of wave function emitted no 
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longer changes with increasing impact parameter. This is shown in Fig.|3] where the fraction 
of wave function emitted Frac = / d 3 r \ip\ 2 is plotted versus the impact parameter. 

For the Au target this percentage starts to saturate around 200 fm whereas for the Ti 
it is around 100 fm and 40 fm for the Be target. At these impact parameters the fraction 
of wave function emitted is around 0.04%. At the saturation, only the noise inherent to 
the numerical method remains and we do not include larger impact parameters in the cross 
section estimate. 

In this figure we clearly see a change in the slope of Frac, between the small and the 
large impact parameters which could indicate the transition between the nuclear and the 
Coulomb perturbation. 




FIG. 3. Fraction of wave function emitted (Frac) versus the impact parameter for the three 
targets (symbols). Dashed lines are fits obtained with the formula air for the large impact param- 
eters (see text). Insert: scaling of Frac — ab@ with A^ 3 (symbols). The solid line is the fit of the 
scaling with e~ lb (see text). 

In order to get a deeper insight in the interplay between nuclear and Coulomb break-up, 
we have studied the mass and charge dependence of the fraction emitted (Frac) with the 
impact parameter. At large impact parameter, above 40 fm, the wave function does not see 
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the target nuclear potential and only the Coulomb field is felt by the core. We observe in 
Fig.|3] that at large impact parameters Frac can be fitted with the formula ab 13 for the Ti 
and Au targets. We found (3 = —2.4 and a which scales quadratically with the charge as 
a = 0.0088 x Z 2 . Fits are shown on Fig.|3| as dashed lines (we assumed the same dependence 
for the Be target). Since the probability of emission is small we can conclude that the 
excitation is perturbative. This is confirmed by the Z 2 dependence of this probability. 

We extrapolated the Coulomb effect to the small impact parameters and by subtracting 
the fit curve to the fraction emitted we obtained a function which we expect to only contain 

2 /3 

the nuclear effect. Indeed, after renormalizing these results with A T , we obtained a quantity 
named Scaling = (Frac — ab l3 )/A^ 3 , which no longer depends on the target (see insert of 
FigfJ). This leads to a cross section that scales with A^ 3 as expected for the nuclear emission 
||17|| . This curve can be fitted in turn by an exponential (Scaling oc c~ lh with 7 = 0.40 ±0.01 
fm- 1 ). 

We would like to emphasize that this empirical fitting method indicates a rather good 
separation between Coulomb and nuclear effects in our calculation. Such a separation, which 
comes as an hypothesis in perturbative framework, is however not natural in our framework 
where both effects are accounted for at the same time and can therefore interfere. By treating 
the Coulomb and nuclear part separately we have controlled that the possible interferences 
are small. 

Using the parameters extracted from the fits and integrating over the impact parameter, 
we extract the variation of the break-up cross section with the charge of the target for both 
the Coulomb and the nuclear interactions. This evolution is presented in Fig.|j. We observe 
a Coulomb cross section that scales with Z im close to the value of 1.725 from ref. |18[ . 



It should be noticed that in the case of the Au target we have checked this analysis by 
turning on and off the nuclear and the Coulomb field separately in the calculation. 



9 



/ Coulomb 



10 Z T 



100 



FIG. 4. Coulomb (squares) and nuclear break-up (triangles) cross sections as a function of the 
target charge (Zt) obtained from the fits to the calculated fraction of the wave function emitted 
and the fit to the Scaling parameter, integrating from the corresponding 6 m , n . The Coulomb 
contribution scales with Z 1,e5 (not shown in the figure). The sum of the Coulomb and the nuclear 
cross-section is also reported as circles joined by a plain line and compared to the experimental 
data (plain circles with error bars) of ref. [3]. 



Once we have the fraction of wave-function emitted we can access the differential cross- 
section by integrating the corresponding probability between b m i n and b max denned in previ- 
ous section. For each impact parameter b n (fin), the cross section is equal to 2nb n Ab n .10~ 2 
(barn) times the modulus squared of the fraction of wave function which is emitted, where 
Ab n is taken to be equal to (6 n +i — & n _i)/2. The step in impact parameter between two 
calculations varied from 1 fm (for small impact parameters) to 40 fm above b = 150 fm. 

The angular distribution of the emitted neutron is extracted by applying the Fourier 
transform to obtain the remnant part of the wave function in the momentum space || 



It should be noticed that to fully take into account the final state interaction this calculation 
should be performed on the asymptotic wave function. We have controlled that p(p) does 



C. Cross sections and angular distributions 



(Fig.0). 




(4) 
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not evolve anymore for longer evolution. Therefore, it can be considered as the final mo- 
mentum distribution of the emitted particles. For small impact parameters, nuclear effects 
are important and the nuclear refraction of the neutron halo gives rise to an emission at 
large angles. This can be seen in the spectra extracted from the calculations performed at 
impact parameters of 10 to 20 fm which exhibit a component at angles above 30 degrees. 
This component is also seen in the density plots of Fig|J and present an anti-correlation 
with the core trajectory compared to the incident direction. This emission to the contin- 
uum is known as Towing Mode. In the case described with stable nuclei, a nucleon from 
the target is pulled out, towed along by the projectile and finally expelled at large angles 
and large velocities in the laboratory frame ||11|| . In the case of the n Be experiments, the 



emitted neutron belongs to the projectile (the n Be) which breaks up as it is perturbed by 
the target nuclear potential. In the n Be frame, this nuclear break-up can be seen as the 
emitted neutron being accelerated by the target potential (see Fig.^). We expect the same 
mechanism to be present for this halo nucleus as for a stable nucleus. However, in the case 
of a halo nucleus, the cross section of this nuclear break-up should be increased compared 
to a stable nucleus, due to the large extension of the neutron wave function and thus the 
larger impact parameter at which it sees the nuclear potential of the other nucleus. This 
will be discussed further when comparing with a more bound wave function (see section D). 
We also display the angular distribution of the initial 2s wave function. In a sudden core 
removal model it would correspond to the distribution of the emitted neutrons. However it 
differs from the presented distributions as well as from the experimental distribution of ref. 
@, showing the importance of the reaction mechanism in the break-up process. 

At larger impact parameters, above 50 fm, only the component between and 15 degrees 
remains (Fig.|5| right). There, the neutron wave function does not feel the nuclear potential 
anymore and the break-up comes from the deviation of the 10 Be core in the Coulomb field. 
Indeed, because the halo neutron of 11 Be is so weakly bound, a light shaking of the core is 
enough to induce the break-up. This, in turn, is not expected with a more strongly bound 
nucleus as it will be discussed further in section D. 
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FIG. 5. Angular distributions of the fraction of the wave function emitted for impact parame- 
ters from 10 to 20 fm (left) and from 50 to 190 fm (right). On the left figure the Fourier transform 
of the initial 2s wave function is also shown as a plain line. 

We have compared our calculation with a time dependent perturbative approach pro- 
posed by Bonaccorso, Brink and Margueron presented in ref. H5|JT9|| (called here perturbative) 
where they treat both the Coulomb and the nuclear break-up as a transfer to the continuum. 
The result we show have been compared with the same optical potential (presented here) 
and the wave function of the neutron halo is taken as the outer part of a 2s wave-function 
and normalized to the wave function used in our dynamical evolution. At small impact 
parameters we know that a large fraction of the wave function is emitted (Fig.[3] shows that 
more than 50% of the wave function is emitted for 10 fm of impact parameter) which in- 
dicates that the perturbative approach may not be adequate in the case of the Au target, 
however at large impact parameters this value becomes small (1% at 30 fm of impact pa- 
rameter) and there the perturbative approach is well justified. We thus expect the results of 
our calculations to be different from the perturbative approach for small impact parameters 
and to become more and more similar as we go at larger impact parameters. Fig.||| shows 
this comparison for three regions of impact parameters for the break-up on a Au target. 



12 




20 40 60 80 



e 

FIG. 6. Angular distribution for the emission of the 11 Be halo neutron using a perturbative 
approach (dashed lines) and our calculation (plain lines) for three regions of impact parameter. 

For the region between 10 and 12 fm the two calculations predict a quite different shape 
of the neutron angular distribution. At large impact parameter, the two calculations ex- 
hibit similar shapes of the neutron angular distributions but a larger cross section is observed 
around 5 degrees in the perturbative approach. This difference, present at all impact param- 
eters, might be due to the remaining differences between the ingredients of the two methods. 
In particular it should be stressed that the wave functions used are not the same in both 
cases. In the perturbative calculation the inner part of the wave function is the Hankel 
function hi while in our approach we use the exact solution of time dependent Schrodinger 
equation. We compared two time dependent theories which show some discrepancies. The 



13 



use of these theories as a spectroscopic tool for halo nuclei requires a better understanding 
of these differences. A deeper analysis is in progress and will be presented in a forthcoming 
article. 

D. Comparison between halo and non-halo neutron emissions 

Our calculation has also been performed to infer the evolution of a strongly bound wave 
function. We used a Wood-Saxon potential with Vo=-70.5 MeV, A=ll, and obtained a 2s 
wave function bound by 7 MeV. The calculation has been performed for a Au target and 
for impact parameters running from 10 fm to 110 fm. The result is shown in Fig. [5] and 
is compared to the break-up of the 11 Be halo for the same impact parameter range. The 
differential cross section of neutron emission for the bound nucleus is more than hundred 
times lower than for n Be below 10 degrees. For large angles, around 40 degrees, the differ- 
ential cross section is about 4 times lower. This latter difference can be understood by the 
extension of the wave function, which is much larger in case of a halo. 
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FIG. 7. Angular distribution for the emission of the 11 Be halo neutron (plain line) and a 
non-halo neutron (dashed line) bound by 7 MeV after the scattering on a Au target summed over 
impact parameters between 10 fm and 110 fm. 
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In particular, it is often argued that Coulomb effects shadow the nuclear mechanism in 
heavy targets like Au. Although the relative proportion of these two effects is largely in favor 
of Coulomb dissociation for a heavy target such as Au, we show that nuclear break-up is of 
major importance for large angle emission with a cross section around 0.5 barn. Whereas 
Coulomb dissociation on a heavy target such as Au could be a direct measurement of the 
binding energy of the particle, the amount of towed particle could in turn bring information 
on the extension of the wave function, and hence answer the question whether we are dealing 
with a halo state or not. Furthermore, in stable nuclei, it has been shown that transfer to the 
continuum due to the Towing Mode might be a tool to infer information on shell structure 



as presented in ref. One might then be able to use this large angle emission to obtain 



additional information on nuclear halo wave-function properties. 

E. Comparison with the n Be data 

To compare our calculation to the data of ref. ||, we extracted the differential cross 
section by summing the calculations from 6 min up to b max . For the Au target we took b max 
= 210 fm, 120 fm for Ti and 40 fm for Be. Calculations are shown in Fig.|, multiplied by 
0.84 to take into account the spectroscopic factor of the 2s state found experimentally and 
reported in ref. [14|. The calculation is compared to the experimental data of ref. || for 



all three targets taking into account the experimental threshold of 26 MeV for the neutron 
detection. Note that our calculation includes all the inelastic channels in which the target 
is excited but also those in which the 10 Be remnant is in an excited state below its neutron 
separation energy. This is also included in the data that measured the 10 Be and the neutron 
in coincidence. However our calculation does not take into account the possible two body 
dissipation since we use a single-particle framework. The calculations (plain lines) are in 
good agreement with the experimental data, both at small and large angles. In Fig.|8| we 
have reported the result of the calculations for three impact parameter regions showing the 
separation between the nuclear and the Coulomb break-up. The simultaneous reproduction 
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of the data for the Au the Ti and the Be targets demonstrates that the Coulomb and the 
nuclear interactions are well taken into account. In these calculations the contribution of a 
neutron in a Id state coupled to a 2+ excitation of the core as deduced from the experiment 



TM has not been taken into account so far. However, since we already reproduce the whole 



cross section for the Au target with 84% of the 2s break-up, the contribution of the Id state 
would seem to be small in this case. 
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FIG. 8. Calculated angular distributions for impact parameter running from b m i n to b max (11 
to 210 fm for the Au target (left figure)) (plain line) for neutron of energy higher than 26 MeV. 
Dots with error bars are data points from ref . [3] . Middle figure is for the Ti target and right figure 
for the 9 Be target. Contribution of the calculation for three impact parameter regions, b m i n to 
15 fm (short dashed lines), 15 to 40 fm (dotted lines) and 40 to b max (long dashed lines) are also 
presented. 



F. Relative energy 

We have calculated the relative energy between the emitted neutron and the 10 Be using 
the momentum of the core after the Coulomb trajectory. This is presented in Fig.|9| for 
three different impact parameter regions for the Ti target. This spectrum shows that large 
relative energies correspond to small impact parameter hence the nuclear break-up whereas 
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large impact parameters lead to small values of the relative energy as expected for the 
Coulomb break-up. 
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FIG. 9. Relative energy between the neutron and the Be core after breakup on a Ti target 
for impact parameter running from b m i n =8 fm to b max =120 fm (plain line). Contribution of the 
calculation for three impact parameter regions, 8 to 15 fm (short dashed lines), 15 to 40 fm (dotted 
lines) and 40 to 120 fm (long dashed lines) are also presented. 



IV. CONCLUSION 

We have investigated the neutron break-up of a halo nucleus in the reactions Au, Ti, 
Be ( 11 Be, 10 Be+n) at 41 MeV per nucleon, in the framework of a time dependent quantum 
model. Results were compared with the experimental neutron angular distributions of Ref. 
and a good agreement was found. Our calculation, that includes both the Coulomb and 
the nuclear interactions, confirms that the forward peaked neutrons are due to the Coulomb 
break-up and that the neutrons emitted at large angles come from the interaction of the 
halo neutron with the target nuclear potential. This diffractive mechanism is called towing 



mode |[n]| . A strong angular correlation between the towed particle and the projectile 
was observed in reactions between stable nuclei and is also expected for the 11 Be break- 
up. However the experiments performed so far did not measure the scattering angle of the 
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remnant 10 Be. Our calculations show that neutrons emitted below 15 degrees arise from the 
Coulomb dissociation and most of the cross section comes from large impact parameters as 
the neutron halo of n Be is weakly bound. The shaking of the 10 Be core by the Coulomb 
field leads to the dissociation of the halo and the emission of the neutron in the forward 
direction. This can also be understood as a Coulomb excitation of 11 Be above the particle 
threshold, followed by neutron emission. This is the only mechanism contributing to the 
break-up when the impact parameter is such that the halo wave function does not overlap 
with the nuclear perturbative potential. The Coulomb break-up is very much hindered for 
strongly bound neutrons, whereas the nuclear break-up decreases by a factor of four due to 
the lesser extension of the wave function. Calculations for the Id wave function should be 
performed for the 11 Be break-up to infer its contribution and understand better the data. 
More generally, those calculations could be used to extract information on the wave function 
of the last bound neutron for unstable nuclei for which the properties are not well known yet, 
provided that a measurement of the neutron angular distribution cross section is performed 
both at large and small angles. 
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